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ABSTRACT 



When low frequency underwater sound interacts with the 
Arctic ice cover, not only will it be reflected from the plane and 
scattered in all directions from roughness elements, but it will 
also be diffracted at leads and reradiated from flexural waves in 
the ice. These phenomena have been studied in an anechoic tank 
by pulse transmission from an underwater point source to a 
series of large floating acrylic plate models, each representing a 
different type of ice cover. The flexural wave speed, the plate 
and lead dimensions and the acoustic roughness are accurately 
scaled, and the specific acoustic impedance contrast is 
approximately modeled by the selection of the acrylic material. 
The physical contributors to the gross reflection coefficient and 
backscattering strength are identified and compared for models of 
a plane ice layer, an Arctic ice pressure ridge, edges of leads, 
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I. INTRODUCTION 



A. PROPAGATION OF SOUND IN THE ARCTIC 

In Arctic waters the existence of a positive sound speed 
gradient in the water column which includes a strongly positive 
near-surface portion which causes a shallow sound duct to be 
formed, forces sound energy to be repeatedly refracted upwards 
to the ice surface. Figure 1.1 illustrates a typical ray diagram 
for sound transmission in the Arctic. Long-range propagation 
measurements in the Arctic Ocean show that the low-frequency 
attenuation and backscattering strengths are more than an 
order of magnitude greater than in the open ocean (Mellen, 

1986) . In order to predict the range at which a noise source 
may be detected in the ice-covered Arctic, it is first necessary 
to understand what happens to the incident sound when it 
interacts with the highly variable ice cover. To do this requires 
a basic understanding of the acoustical characteristics of the ice 
cover . 

B. DESCRIPTION OF THE ARCTIC ICE COVER 

In general, the Arctic ice surface is variable in thickness, 

density, elastic constants and shape. As illustrated in Figure 
1.2, the sea ice can be partitioned into the following acoustically 
significant categories: level ice (70 %) , pressure ridges and 
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Figure 1.1. A typical sound velocity profile and associated ray 
diagram for underwater sound propagation in the Arctic 
(Diachok, 1974) . 



INTERACT ICH OF LOW FREQUENCY SQUWP WITH JH£ ARCTIC ICE CAEQei 
■£QTE?tIlAL C flilTRlSUTORS 



AIR 



ICE SAIL 



FLEXURAL WAVES OW ICE 




0 s Diffraction from leading end trailing edges of plates 
K = Keel scatter from deterministic ridges (ice keels) 

S r Scatter from steep-sloped, random! y-rougn, rut>oie fields 



B = Boundary wave at sieeo-sloped, randomly- rough, ruOOle fields, near grazing forward scatter 
Suosonlc 



F = Flexural (surtece) wave, SuOsomc, evanescent in water, radiating into the airsoace 
aoove the plate 



Figure 1.2. Scattering of low frequency sound from the 
Arctic ice canopy. Acoustically significant features of 
the ice canopy include level ice, pressure ridge keels, 
hummock fields and edges of leads and polynyas. 
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rubble fields (25 %) , and leads (5 %) . The number in 
parentheses is the approximate percentage of the Arctic ice 
cover each category represents (Weeks et al, 1971). This 
experiment is concerned only with kh < 1, where h is the ice 
thickness and k is the horizontal wavenumber. Each of these 
ice categories is briefly described below. 

1. Level Ice 

Sea ice is termed level if its surface changes less than 
1 part in 40 (Wadhams et al, 1985, p.1072). The thickness of 
level ice depends on the age of the ice, with first year ice 
generally less that 2 m thick and multi-year ice having a mean 
ice thickness of 2.5 to 3 m (Ackley et al., 1976, p. 9). 

2. Pressure Ridges 

Ridges are formed by the interaction of adjacent ice 
floes. The type and size of a ridge is a function of the thickness 
of the interacting ice and the degree of compressive or shear 
forces being exerted on the ice floes. 

Hibler et al (1972) compiled above-ice data from visual 
and laser profilometry and underice data from submarine 
underice profiles to develop statistical models of pressure ridge 
keel depths and the spacings between ridges. To insure that 
only pressure ridges are dealt with in the above ice and 
underice profiles, a cutoff height (or depth) was used, i.e., 
ridge height and keel depth distribution data are reported from 
a set that includes all ridges greater than a given cutoff. The 
average pressure ridge keel depth reported by Hibler et al 
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(1972) from the set of all ridges greater than 6.1 m, was 
9.6 m in the Central Arctic and 11.4 m in the Canadian 
Archipeligo. They further report that the ratio of sail height to 
keel depth is on the order of 1:5. They report the average 
spatial distribution of pressure ridges as 4 . 3 ridges/km in the 
Central Arctic and 9.6 ridges/km in the Canadian Archipeligo. 

Hibler et al (1974) reported that ridging intensity is 
seasonally dependent. They determined that the spatial 
distribution of ridges in the Beaufort and Chukchi Sea regions 
was 1.83 ± 0.6 ridges/km in March, and 0.99 ± 0.61 ridges/km 
in October and November. Furthermore, they observed that 
all the ridges investigated had similar slope angles, such that 
the ratio of the base width to the ridge height was 
approximately four. 

A specific multi-year pressure ridge was the object of 
investigation by Kovacs et al (1973) . This particular ridge had 
a sail height of 4 m and a keel depth of 13 m. The multi-year 
pressure ridge differs from the first year pressure ridge in that 
the latter is poorly bonded and contains many spaces between 
the blocks of ice that comprise the ridge. In a multi-year ridge 
these voids eventually are filled with ice due to the seasonal 
melting and re-freezing of the ice. 

Wright et al (1978) report values for the keel depth 
to sail height ratio and keel width to keel depth ratios slightly 
different from those above, i.e., 3.2:1 for keel to sail ratio and 
3.3:1 for the width to keel depth ratio. They also report that 
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keel depth statistics show little or no seasonal variability 
probably due to the fact that the underice surface does not feel 
the effects of the ablation of the ice surface during the warming 
months. 

In the present model the ice ridge is assumed to be a 
triangular cross sectional (prismatic) wedge below the plate and 
a similar smaller prism above the plate. The fine structure is 
omitted in this low frequency model. 

3. Hummock Fields 

In areas of severe shear forces ridges can be pushed 
together to form a series of adjacent ridges which would appear 
from above like a farmer’s plowed field (Hibler and Ackley, 

1975) . This type of non-random rubble field occurs where the 
forces acting on the ice are directional such as in a narrow 
strait or butted against a land mass. In another situation 
where the forces are less constrictive, such as in the Central 
Arctic Basin the ridges could be pushed together to form a 
nearly random field of jumbled blocks interspersed with large, 
nearly vertical ridge segments. It is this latter description 
which was modelled for this experiment. 

4. Leads and Polvnas 

An open lead or polynya in the ice cover is an area 
where there is no ice. A covered lead exists where the 
thickness of the ice is much less than the that of the 
surrounding ice floes, i.e., below 1 m thickness 
(Wadhams and Horne, 1980). Leads are generally considered to 
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be cracks in the ice where the length of the crack is much 
greater than the width. Polynyas on the other hand are large 
open areas up to several kilometers wide. For this experiment, 
open leads of 6, 12 and 18 m width have been modeled. A 
single edge of an open polynya has been modeled by using a 
single plate edge. Chapter III describes these models in greater 
detail. 
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II. RESEARCH FACILITIES AND EQUIPMENT 



A OCEAN PHYSICS LABORATORY AND ANECHOIC TANK 
The anechoic tank provided a low noise environment 
sufficient to conduct this experiment. The dimensions of the 
anechoic tank are: length - 7.3 m, width - 1.8 m, depth - 
2.2 m Aluminum-impregnated rubber conical sound absorbers 
line each wall. Pulsed signals were transmitted and the signal 
return was identified, isolated and analyzed utilizing a judicious 
choice of a sampling window at the signal processing end. 

Figure 2.1 illustrates the equipment set-up utilized throughout 
the experiment. An equipment list is provided in Table I. All 
equipment used is available on the commercial market with one 
exception - the Random Surface Measuring Device (RSM) . This 
device, which was designed by a former NPS student LT. S. 
Kasputis (1984), allows the accumulation of a large data base 
for statistics on a given surface. It can determine heights on 
slopes of less than 70 degrees with a precision of 20 pm. 
Additional equipments which merit a more detailed description 
are discussed below. 



19 



HP-467 A 
Pvr Amp 




= 1 MHz 



Figure 2.1: Schematic diagram of the signal processing 

flow path. 
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B TRANSMIT AND RECEIVE SYSTEMS 

1. Hewlett Packard Function Generator, Model 3314A 

The 3314A Function Generator is a digital, multi-mode, 
programmable function generator which features preset Sine, 
Triangle and Squarewave functions and programmable arbitrary 
waveforms from ImHz to 19.999MHz and a maximum output 
voltage of 10 volts peak-to-peak. The function generator was 
utilized in the N CYCLE MODE which outputs a burst of N 
complete cycles, starting when an internal or external trigger 
edge is received. The excellent stability of the function- 
generator permitted repeated data runs in minimal time. 

2. Nicolet FFT Wave Analyzer . Model 660B 

The Nicolet wave analyzer is programmable and performs 
its own analog to digital conversion. The sampling rate is 
limited to 250,000 samples pen second. It has a built-in CRT 
screen that allows real time observation of either the received 
time signal or its frequency spectrum. This equipment was 
utilized in parallel to the signal processing equipment to permit 
initial detection of the signal of interest. 

3. Interface Technology Timing Simulator. Model RS-648 

The timing simulator was used to trigger the function 

generator, the Nicolet and the signal processing equipment. It 
has a 50 nanosecond timing resolution. 
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4. Celesco LC-10 Hydrophone 

LC-10 hydrophones were used as both transmitter and 
receiver. The specifications for the experiment required a point 
source (ka « 1, where 'k' is the wave number and 'a' is the 
linear dimension of the transducer) and receiver with rise and 
decay times sufficiently short to permit the reception of 
undistorted pulsed signals in an environment where the time 
window was critical. The time duration of the pulse is 
important since all of the desired signal must be sampled in a 
time window free of the unwanted arrivals from the tank or 
plate reflections and diffractions from the edges of the model. 
Additional requirements were a relatively high and flat 
free-field voltage sensitivity and an omnidirectional directivity 
pattern in the frequency range of 30 kHz to 100 kHz. The 
LC-10's directivity pattern was omni-directional to within ± 2 
dB. The greatest disadvantage of the LC-10 as a source is its 
low output power, however, the combination of short ranges 
(less that 2 m) and the low ambient noise level in the anechoic 
tank permitted its use with excellent results. 

C. SIGNAL PROCESSING EQUIPMENT 

Signal processing was accomplished using an IBM PC/XT 
personal computer fitted with a Computerscope Model ISC- 16 
Signal Processing system. The Computerscope ISC-16 system is 
a fully integrated hardware and software package designed to 
permit the IBM to perform as a data acquisition and analysis 
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laboratory instrument. The ISC- 16 package consists of a 
16-channel analog to digital converter board, an external 
Instrument Interface and Scope Driver software. The package is 
capable of receiving 16 channels of data input at an aggregate 
sampling rate of up to 1MHz. Digital conversion is achieved 
with 12 bit accuracy over an input range of -10 to +10 volts 
The Scope Driver software permits the IBM CRT to emulate a 
digital storage oscilloscope with numerous waveform 
manipulation and display modes, and to output high resolution 
graphics to a dot matrix printer or X/Y plotter. Figure 2.2 is a 
graphic produced by the Computerscope system and illustrates 
the systems excellent signal processing capabilities. The package 
is capable of directly accessing 640K of RAM. 
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TABLE I 



EQUIPMENT LIST 



Abbreviation 


Nomenclature 


Timing Simulator 


Interface Technology Timing 
Simulator/Word Generator, Model 
RS-648 


Function Generator 


Hewlett Packard Function Generator, 
Model HP-3314A 


Power Amplifier 


Hewlett Packard Power Amplifier, 
Model HP-467A 


Source or Receiver 


Celesco LC-10 Hydrophone 


Pre-Amplifier 


Ithaco Low Noise Pre-Ampifier, 
Model 1201 


Nicolet 


Nicolet Dual Channel FFT Analyzer, 
Model 660B 


Voltage Amplifier 


Hewlett Packard Voltage Amplifier, 
Model HP-465A 


Signal Processor 


IBM PC/XT with Computerscope 
ISC-16 Data Acquisition and Analysis 
Package 


Scope 


Kikusui Oscilloscope, Model COS 5060 
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Figure 2.2: 

Time domain voltage response to an input sinusoidal 
pulse, driven at 7.07 volts RMS. Sinusoid frequency 
is 56 kHz. 



Input pulse: y(t) = A sin( 2 tt ft ) [1[t) - 1(t-2/f)] 



Peak Identification: 

QQ - Direct signal 

(2) — Surface reflected signal 

Q>) - Specular scatter from sides of tank 

( 4 ) - Diffraction from corner of watei — air interface 

and side of tank 

( 5 ) — Specular scatter from the bottom of tank. 
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III. DESCRIPTION OF PHYSICAL MODELS 



Before construction of the "ice" models could begin, it was 
necessary to select the material which best approximated the 
physical characteristics of Arctic ice. The internal structure of 
sea ice, such as its saline content and the distribution of melt 
inclusions, is beyond the scope of our interest; furthermore for 
the low frequency regime which is our interest, average values 
provide a satisfactory acoustical description. Therefore, for this 
experiment it was necessary to emulate the bulk properties of 
the ice cover. Acrylic plate was chosen because it has several 
important physical properties that parallel closely those of Arctic 
sea ice. Specifically, the longitudinal wave speed and flexural 
wave speed (for scaled thicknesses) , Poisson Ratio, and the 
characteristic impedance (pc) for acrylic plates all fall well 
within the measured range for Arctic sea ice. The density of 
the acrylic plate is about 25% too high and the shear and 
compressional attenuations of the plates are somewhat higher as 
well. Table II summarizes the comparison between Arctic ice 
conditions and our laboratory acrylic "ice sheet" floating in a 
fresh water tank. 
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TABLE II 

COMPARISON OF THE PHYSICAL PROPERTIES OF 
ARCTIC ICE TO THE ACRYLIC MODEL LABORATORY 
CONDITIONS 



PROPERTIES 


ARCTIC 


LAB 


MODEL/ICE 










RATIO 


CHARACTERISTIC 


RANGE TYPICAL 


MODEL 




1. YOUNG'S MOD 


1. 7-9.0 


4.1 


4.1 


0.46-2.41 


E, (Pa x 10 9 ) 


(1) 


(2) 


(3) 




2. DENSITY 


760-960 


930 


1148 


1.20-1.51 


P2 (kg/m 3 ) 


(4) 


(4) 


(3) 




3. POISSON'S RATIO 


0.30-0.36 


0.33 


0.316 


0.88-1.05 




(5) 


(5) 


(3) 




4. BULK COMPRES- 


2800-3500 


3150 


2240 


0.50 - 1.45 


SIONAL WAVE 
SPEED (m/s) 


(6) 


(6) 


(3) 




5. COMP. WAVE 


0.3F-0.05F 


0.06F 


0.424F 


1.41 - 8.48 


ATTENUATION {a) 
(dB /m), F(kHz) 


(20) 


(20) 


(3) 




6 BULK SHEAR 


1500-1860 


1700 


1373 


0.74 - 0.92 


WAVE SPEED (m/s) 


(6) 


(6) 


(3) 




7. SHEAR WAVE 


1.8F-0.30F 


0.36F 


0.853F 


0.47 - 2.37 


ATTENUATION (o) 
(dB/m), F(kHz) 


(20) 


(20) 


(3) 
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TABLE II (CONT'D) 

COMPARISON OF THE PHYSICAL PROPERTIES OF 
ARCTIC ICE TO THE ACRYLIC MODEL LABORATORY 
CONDITIONS 



8. BULK CHAR. IMP. 1.06-3.54 

P2 C 2> ( 8 ) 

(kg/m 2 -s) X 10 6 

9. FLEXURAL WAVE 252-2357 

SPEED (m/s) (11) 

10. FLEXURAL WAVE Data not 

ATTENUATION available 

(dB/m), F(kHz) 

11. PLATE THICKNESS 0. 5-6.0 

h (m) (4,7) 

12. THICK/WAVELENGTH 0.007-0.1 

h/X (14) 

13. VEL. RATIO 0.97-2.48 

C 2 /Ci (17) 

(ICE/WATER) 



2.07 2.57 0.73-2.42 

(9) (10) 



898-1270 


700-1150 . 0.49-2.78 


(12) 


(3) 




1.5F 




(3) 



2.0 


0.00318 


1:600 


(7) 


(3) 




0.09 


0.04-0.11 


0. 4-1.2 


(15) 


(16) 




1.54 


1.51 


0.61-1.56 


(18) 


(19) 
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